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TECHNICAL MEMORANDUM X-G4843 


A PROCEDURE FOR CALCULATION OF BOUNDARY LAYER 
TRIP PROTUBERANCES IN OVEREXPANDED ROCKET NOZZLES 

SUMMARY 


Side forces occur during rocket engine start and shutdown in ambient 
pressures greater than zero. These forces are caused by non-axlsymmetrical 
flow separation. One method of avoiding this undesirable efitect Is by using 
a boundary layer trip to force a clean, axisymmetrlcal separation line. 

Since boundary layer disturbances generate shock waves and reduce 
nozzle performance, boundary layer tilp sizing must represent a compromise 
between performance reduction and trip effectiveness. Procedures for 
calculating trip sizes, required positions, scaling effects and performance 
reduction plus a comparison between theory and available experimental data 
are presented. 


INTRODUCTION 


During rocket engine start and shutdown in ambient pressures greater 
than zero, nonaxlai thrust components act as side loads on rocket nozzle 
walls. These unsteac^, randomly oriented side forces are caused by un- 
symmetrlcal flow separation, and can impose a requirement for additional 
structural support thereby Increasing the weight of the engine and mountings. 
Since such weight increases are undesirable, several methods of reducing the 
magnitude of the side forces have been proposed. [ 1] 

One method of reducing side loads is to force clean, axisymmetric 
separation lines. This can potentially be accomplished through the us I 
boundary layer trips such as secondary injection, wall angle discontinuk;.ies. 
circumferential grooves, and circumferential trip protuberances. The use 
of a protuberance can be accomplished by installing a trip wire which has an 
advantage that it can usually be installed after the development of an engine 
and nozzle without any change of the basic configuration. 



Any protuberance of a nozzle wall generates drag losses and shock waves 
which reduce engine performance. Since the size of the protuberance or trip 
wire affects the flow separation behavior, performance reduction and flow 
separation characteristics are related. A trip wire sizing and positioning 
procedure Is required, therefore, to provide the combination of precise 
separation and minimum performance reduction. 

This report describes one method for sizing, scaling, and positioning 
of trip protuberances. Some experimental results are presented, allowmg a 
comparison between theory and experiment. The main purpose howev^,r. Is 
to outline the rationale of the proposed calculation method. 


TRIPWIRE SIZING 


The objective In using trip wires Is to achieve a forced separation of the 
flow. This separation should always occur at a higher wall pressure than 
without a trip wire. One or more trip wires in a nozzle are used as shown 
In Figure 1. During transient start and shutdown, the chamber pressure Is 
lower than during main stage operation and the flow separates from the wall. 
At a given chamber pressure, the gases expand into the nozzle and flow over 
the trip wires, 1 to I - 1. At trip wire I, the flow disturbance is sufficient 
that the flow separates. Without a wire Installed, the flow would have sep- 
arated further downstream In the nozzle. When the chamber pressure is 
Increased, the gases finally overflow this trip wire. Then the next wire, l+l, 
has been positioned such that at l+l forced separation is achieved. A further 
Increase of the chamber pressure lets the flow jump to wire i+2 and at last the 
nozzle flows full. 

Although this approach of having many or an almost "continuous" series 
of trip wires along the nozzle wall seerrs to be promising, some of the disad- 
vantages accrued appear to favor the use of oniy one or a few trip wires. The 
performance loss Increases with the number of trip wires, and the change of 
the flow field over the trip is so severe that a certain distance between the 
trip wires is required so that the disturbances can decay. 

The flow separation process is a boundary layer effect. Therefore, a 
short description of the boundary layer model, which Is used for calculations. 
Is necessary. 




Figure 1. Nozzle flow separation, with and without boundary layer trips. 


Boundary Layer Model 

The boundary layer which approaches the trip wire is presented in 
Figure 2, together with the associated nomenclature. Since performance 
efficiency requires a small trip size, a trip height of about 1/10 of the boundary 
layer thickness should not be greatly exceeded. Therefore the distribution of 
flow properties at the boundary layer ec^e is of little interest. 

The usual assumption for the velocity profile within the boundary layer 
is a 1/7 power law. Denoting the velocity as u, the velocity at the boundary 
layer edge as u , the boundary layer thickness as 6, and the coordinate normal 

to the wall as y, this distribution is e^qpressed by: 
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Figure 2. Boundary layer model. 



( 1 ) 


The existence of a laminar sublayer will be neglected within the accuracy of 
this calculation. 

The temperature distribution can be related to the velocity distribution 
by the Crocco equation: 



( 2 ) 


where, T^^ is the local stagnation temperature, T^ is the stagnation tempera- 
ture at the boundary layer et^e, and T is the wall temperature. Since T is 

w r 

only slightly different from the combustion chamber temperature T , with the 
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energy equation 
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the local static temperature T can be obtained from: 


T = 



T \+ T 

w 


u‘ 



( 4 ) 


Rearranging with the static temperature at the boundary layer edge 



ib) 


yields 



In these equations, M describes the Mach number at the boundary 

layer edge and y the isentropic exponent of the gases, which are assumed to be 
ideal. 


Upstream of the trip outside of its influential range a constant static 
pressure across the boundary layer is assumed. With the gas equation the 
density p is expressed by 



( 7 ) 


where p represents the static pressure and R the gas constant. 



Performance Loss 


The estimation of the performance loss by the trip protuberances 
requires some assumptions. The trip wire is small ompared with the bound- 
ary layer and therefore the shock, which is generated in front of the disturb- 
ance, becomes weaker towards the boundary layer edge and should disappear. 
Then the wall pressure further downstream of the trip deviates only slightly 
from the undisturbed wall pressure and the performance loss due to wall 
pressure reduction can be neglected. With this assumption, the performance 
loss is generated only by the dri^ of the trip wire. 

The drag is normally calculated by the dynamic pressure and a suitably 
defined drag coefficient C^. Since the velocity and density upstream of the 

trip varies normal to the wall, an integration is necessary. Ihen the drag per 
length is (Fig. 3) 


tw 

Df = / . (8) 


In all effects such as low pressure at the back side of the trip protiiberance 
are included. With the equations (6) and ('') , the drag can be expressed by 



In order to simplify the calculations^ a constant, properly defined 

is useful. An approximate value can be obtained from the various drag meas- 
urements of different body shapes. Rearranging with the boundary layer ed^^c 
velocity of sound and equations (l) and '' ' results in 
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( 10 ) 


The exact integration of equation (lO) requires the use of a decompc.->ition to 
partial fractions. An approximate solution can be obtained by the fact that 
above a distance of 0.01 6, the integrand of (lO) increases rather linearly. 

NOZZLE AXIS 



Therelore, equation (lO) can be written with little error 



The factor k represents the rather linear increase of pu with y/6 in the region 
discussed and has a value of 0.5. 

For a quick estimation of the second term on the right side of equation 
(ll) , some typical numbers for the different symbols are; 

M = 3 
e 

d ^ /6 = 0. 1 

tw 

T /T = 0.5 
w c 

y = 1.25 
and equation (ll) yields 


0.35 pyd 

n = n M" ^ 


( 12 ) 


The thmst loss is obtained by multiplying the drag per unit length with 
the nozzle perimeter, using only the axU'l component. With the nozzle diam- 
eter d^ and the wall angi. 0^, the thrust loss can be expressed by 

AF. = D, ird cos ^ . (13) 

tw 1 n n ' ' 












The previous listed typical values can be used for a simple, approximate 
equation of the thrust loss. Assumir.g a drag coefficient of 0.8, which is sim- 
ilar to those for spheres or cylinders, equations (l2) and (13) result in 


AF_ = 0.5 p d^ d 
tw e tw n 


(14) 


Equation (14) indicates that only very small trip sizes lead to negligible per- 
formance loss. 


Wire Size and Flow Separation 

A turbulent boundary layer in a rocket nozzle can only withstand a 

certain overexpanded condition. This condition is expressed by the separation 

criterion, which describes the ratio of the minimum full flowing wall pressure 

p. to the ambient pressure p . The use of a trip :;hai^es this behavior, always 
1 8 

achieving separation at a higher wall pressure than would be achieved without 
a trip. 


Although the natural separation phenomenon in an overexpanded rocket 
nozzle is far from being solved completely, cer+^in methods are available 
which describe the experimentally observed res s reasonably well. * The 
theory', which was derived by Crocco and Probste.n (2l , achieves the best 
agreement with test data. The numerical results of this theory and the avai- 
lable separation data of hot firing nozzles are plotted in Figure 4. 

The calculation of the relation between trip protuberance size and 
separation behavior can be done by a momentum approach.^ The control 
volume, which is used tor the analysis, agrees with that of Crocco-Probstein 


1. Schmucker, R. H. : Status of Flow Separation Prediction in Liquid 
Propellant Rocket Nozzles. NASA TM X (to be published) , Marshall 
Space Flight Center, Alabama 35812. 

2. Another method for the determination of size effect on separation behavior 
is presented in the Appendix. This approach is based on the suggestion that 
the shock waves in the boundary layer change within a small region of the 
wall pressui*e, so that separation can occur. 
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Figure 4. Experimental flow separation data and Crocco-Probstein flow sepa- 
ration criterion (Data points represent averaged hot firing data) . 

(2] , and is presented in Figure 5. The momentum balance between the point i. 
which is upstream of the trip, and a point downstream of the trip (which is 
described by the subscript a for convenience) , yields: 



where I is the momentum and the subscript wt describes the condition ”with 
trip". Without a trip, equation (15) is written 



( 16 ) 
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Figure 5. Control volume for trip size calculation. 

with subscript ot as ’’without trip”. The latter condition corresponds to 
normal separation in a nozzle. Combining equations (15) and (16) results in 





(17) 


The pressure in equation (ll) for the protuberance drag is the un- 
disturbed pressure upstream of the trip. Relating the wire drag in equation 
(17) to this pressure results in 
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The separation conditions with a trip protuberance are calculated by using 
equation (18) , the -boundary layer data, and the separation data without trip. 
The latter can be obtained from Figure 4. 

Equation (l8) indicates a rather linear relation between the separation 
pressure change p. - p. and the wire size, normali::,ed with the boundary 
wt ot 

layer thickness. With increased wall cooling, the effectiveness of the trip 
increases since the density of combustion gases near to the wall is increased. 

The combination of equation (18) and (l2) , together with numerical 
values allows the derivation of a simple, approximate equation for a rough 
estimation of the trip effect 



M. is the Mach number at the boundary layer edge at the trip location (at point 


The inversion of equation (18) allows the calculation of the trip size for 
a desired separation pressure change. Since the trip drag of equation (ll) 
depends on the 1/7 and 2/7 power of a trip size, only an iterative solution is 
possible. A simple approximation is obtained by equation (19) , using lor 
the pressure increase. 



Equation (20) shows a strong influence of the Mach number, indicating a 
decrease of the normalized wire size along the nozzle wall. But since the 
boundary layer thickness increases with the lei^th, the absolute value of the 
trip size increases also. 


! 


In Figure 6, some experimental data together with theoretical calculation 
are presented. The test data are obtained from a 4k H^/LOX engine at NASA's 
Marshall Space Flight Center. The}' represent the minimum wall pressure 
p. during the test. Since this nozzle was uncooled and therefore changed its 
wt 

temperature during the tests, the calculations were performed for two wall 
temperature ratios. 


0.146 0.30 





Figure 6. Calculated and measured separation pressure ratio (ML=3.6,7=1.26) 
(Test data obtained with NASA-MSFC 4-k H 2 /O 2 engine) . 

The agreement between theory and experiment is rather good. The 
theoretical calculated values represent an upper limit, which is not exceeded 
by the test data. The measured separation pressures are normally lower than 
the theoretical limit, since the trip location does not agree with the upper trip 
position. (See later discussion of Trip Protuberance Position.) The test data 
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with the smaller trip size asree with the theoretical line. This indicates that 
the test results should be very sensitive to boundary layer changes. During 
these tests the separation point jumped from one trip to the next, and was 
probably caused by the wall temperature increase. 


SCALING OF TRIP PROTUBERANCE 


The performance loss and effectiveness depends on the drag of the trip 
protuberance. The exact numbers and relations of drag coefficient, velocity, 
and temperature distribution are not known. The measurement of separation 
behavior avoids this difficulty, since integral quantities in equations (14) and 
( 20 ) are obtained. 

For a measured separation pressure ratio and a fixed Mach number and 
trip shape, equation (20) can be simplified and the trip size is 





(21) 


The trip protuberance height is a certain fraction of the boundary 
layer thickness. If the boundary layer thickness is known, the required size 
is obtained from equation (21) using test results for Ci- 

If no boundary layer data are available, a rough estimation is possible 
using the boundary layer thickness of a flat plate. For convenience a conical 
nozzle is assumed (Fig. 7) . Then, with the length x which starts at the throat, 
the boundary layer thickness can be written [3l 


0.2 



where is the velocity at the boundary layer edge, 5. is the boundary layer 

thickness at the trip location i, and v is the viscosity. The length x can be 
approximated by 
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(23) 


n 

2 sin 0 
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Introducing the thrust equation 


F 


d2- 
t 4 



(24) 


where d describes the throat diameter and C the thrust coefficient, thus 
t F 

yielding; 


/ p \ 0 .5 g 0 .5 

VC/ C_°‘®sin e IT®’® 
F n 


(25) 


Combining (22) and (25) results in 



t 


(26) 


where all constants are comprised in C 2 . The trip wire size increases with 
the 0.4th power of thrust and decreases in the same way with chamber 
pressure. 

The previous described scaling procedure is cnly valid for scaling to a 
different engine with the same boundary layer conditions. A scaling to other 
conditions, such as expansion ratio c, Mach number, etc., requires a change 
of the constant C 2 * These effects can be estimated by the proportionality 


(27) 
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NOZZLE AXIS 



Figure 7. Geometry of a conical nozzle. 

but the accuracy of such an approach is very questionable. Equation (27) 
should therefore only be used for the estimation of trends. 

TRIP PROTUBERANCE POSITION 


The calculation of the trip protuberance size requires, besides the 
wall temperature and the Mach number, the desired wall pressure at which the 
flow separates. It is obvious that this condition can only be fulfilled at a cer- 
tain chamber pressure. During transient conditions, the chamber pressure 
changes and the trip device is only effective during a certain chamber pressure 
range. Therefore, the trips have to be positioned such that natural separation 
cannot occur anywhere. 

For calculation of the possible location range for the trip protuberance. 
Figure 8 can be used. At a certain chamber pressure, an undisturbed wall 
pressure profile is obtained like that of Figure 8. According to the wall pres- 
sure distribution, a Mach number can be calculated at every station and, with 
the data of Figure 4, the natural separation pressure can be plotted. The inter- 
section of the wall pressure profile and the separation pressure profile de- 
scribes the natural separation point. Since a deviation from these ideal values 
always occurs, a small scatter must be introduced to get reliable numbers. 
Assuming a certain trip size allows, together with the local boundary layer 
thickness and wall temperature, the calculation of the trip drag. With the 
separation pressure and equation (18) , this leads to the separation pressure 


with trip. The intersection of this distribution with the undisturbed wall 
pressure results in the upper limit for the trip position. Placing the trip up- 
stream of this point would result in a reattachment of the flow after the trip. 
The use of a trip protuberance downstream of the lower limit does not effect 
the natural flow separation at all. Therefore, a trip wire of a fixed size 
should only be placed in this described range. 



POSITION 

AXIAL DISTANCE PROM THROAT 


Figure 8. Position of trip protuberance (d = const) . 
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A changing chamber pressure varies only the undisturbed wall pressure 
profile. If the trip is positioned in the previous mentioned region of the 
nozzle, an increasing chamber pressure moves the condition to the upper 
limit. Finally the chamber pressure is high enough that the gases overflow 
the trip. In this case, at the lower limit of the trip location, a new' trip de- 
vice has to be located so that natural separation can be avoided. This holds 
true for every trip protuberance. The nozzle exit can also be considered as a 
trip. 

The separation pressure increase with trip depends on the trip size. 
Smaller trips are less effective and require a larger number of devices. 

Bigger trips allow larger spacing and are not so sensitive to boundary layer 
changes. This favors the suggestion that for flow separation with trips, only 
a few or even one trip wire should bo used. 


CONCLUS ION 

Boundary layer trip protuberances are a promising approach for side 
load reduction in overexpanded rocket nozzles by forcing a clean, axisym- 
metrical separation line. Effectiveness can easily be verified by installing a 
circumferential trip wire. 

A sizing, scaling, and positioning procedure for the trip protuberances 
is described, and is based on the momentum equation. This derivatiop 
indicates that performance loss and trip effectiveness are related. The 
performance loss is very sensitive to trip size, but only very small trips are 
necessary for the significant change of the separation behavior. Some ex- 
A perimental data are presented, which show good agreement with the theoretical 

; results. This suggests that the theory presented describes the real phenomena 

rather well. 

The scaling of trips requires some information about the boundary 
f layer. A simplified scaling procedure relates thrust, chamber pressure, and 

I trip size for a fixed boundary layer condition. Positioning and spacing of trips 

i : is based on the natural separation behavior and the change caused by the tri))s. 

j Larger spacing, which is desired for less boundary layer disturbance, can be 

^ : achieved by a few biggp'* trips. The last trip has to be in such a position that 

I ' the nozzle exit is being considered as a final trip. 
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-PEND IX 

CAirUi ^ ;>l OF THE TRIP SIZE BY 

WALL . iJRE CHANGE DUE TO SHOCKS 


The ilr,w fiel . ar oui-.d the trip protuberance and the wall pressure 
distribution is presented in, Figure Al. Upstream of the wire the flow is 
compressed and si-t'.i rates from the wall. A curve shock emerges in this 
region, decreasing in strength with Increasing distance from the wall. The 
separated stream line reattaches at a certain point from the protuberance. In 
the separated region the wall pressure increases. Downstream of the first re- 
attachment point the flow expands and the wall pressure drops. The trip size 
is so small that the area ratio of the nozzle can be considered as unchanged. 

Due to the shock wave the local stagnation pressure in the boundary layer 
drops. At a certain point from the trip protuberance, the flow separates again 
and the wall pressure drops. If the chamber pressure is high enough, the 
gases flow over this protuberance and reattach downstream of the trip. There- 
fore, a reattachment shock is generated, which leads to a further stagnation 
pressure loss. Since the disturbance of the flow field occurs only in a small 
region near to the wall, the mixing process damps the differences of the 
stagnation pressures and the deviation from the theoretical wall pressure 
profile decays very fast. The final pressure distribution should be slightly 
below the disturbed pressure profile, but the deviations are almost too small 
to be measured. (Wall pressure measurements in a small nozzle w'ith rectan- 
gular trip wires do not show any deviation from the wall pressure without 
trips.) If the chamber pressure is low enough, the expansion downstream of 
the trip prevents a reattachment and a forced separation is achieved. 

Downstream of each reattachment point the static wall pressure is 
lower than upstream of the wire since shock waves decrease the local stag- 
nation pressure. This favors the suggestion that with a fixed trip configuration, 
a certain chamber pressure range exists which will always lead to forced 
separation. The lower limit of this range corresponds to separation at the 
edge of the trip.^ The higher limit corresponds to the condition just before re- 
attachment after the trip. 


3. The flow always separates at the trip or upstream of the trip if the 

chamber pressure is lower than this lower limit. Therefore this chamber 
pressure must be considered as a theoretical lower limit. 
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WAU. PRCSSURE 


prcscurc change is described by a normal shock. This can be considered as 
a mean between the previous mentioned upper and lower limit of the separation 
condition. 



Figure A2. Simplified pressure distribution at trip. 

Denoting a,s the stagnation pressure ratio of the trip edge stream 
lino across the shock waves, the expression for the wall pressure is 

f 
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The stagnation pressure ratio is obtained by the usual shock relations (4) , 


Y 



In the interesting Mach number range this expression can be simplified by a 
linear relation. 


= 1.45 .. M. 0.36 . ' (A3) 


The Mach number is calculated by equation (l) , (6) , and (7) t 



(A5) 


M. represents the boundary layer edge Mach number at the trip position. 

A comparison of the results ot equation (A5) with the test data of 
Figure 6 indicates that the assumption of a normal shock leads to discrepancies 
witn the measurements. Althougn equation (A5) predicts rather small trip 
«?i 7 .es and the correct trend of the wall temperature effect, the deviation from 
the real phenomenon is excessive. Therefore trip wire sizing based on shock 
generation is not an acceptable approach. 
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